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O . ABSTRACT 
<N 

<j_) , The apparently repeating microlensing event OGLE-2003-BLG-095 is ana- 

lyzed. Data were obtained from the OGLE Internet archive and exist in the 
public domain. The source is relatively bright, with an unmagnified (but possi- 
bly blended) J-band magnitude of 15.58, and the signal-to-noise ratio of the data 
is excellent. The light curve shows two distinct, smooth peaks characteristic of 
a double microlensing event. It can be modeled as either (1) microlensing by a 
^ ! binary lens or (2) microlensing of a binary source, with the latter model providing 

a statistically superior fit. However due to apparent low-amplitude variability of 
the source, the interpretation is somewhat ambiguous. OGLE-2003-BLG-095 is 

O 

2 ' Subject headings: binaries: general — gravitational lensing 



only the second possible case in the literature for microlensing of a well-resolved 
binary source. 



1. Introduction 

Gravitational microlensing surveys were originally suggested by Paczyhski (1986) as 
a means of detecting dark matter in the Galaxy in the form of massive compact objects, 
commonly abbreviated as MACHOs following Griest (1991). Various efforts were undertaken 
to search for such objects, including the Optical Gravitational Lensing Experiment (OGLE; 
Udalski et al. 1992), the MACHO project (e.g., Alcock et al. 1993), and several others. 
Currently ongoing projects, most notably OGLE and MOA (Bond et al. 2001), are detecting 
hundreds and dozens of new microlensing events each year, respectively. The new events 
found by these groups are made public in real time to maximize scientific gain by enabling 
follow-up of interesting objects by the astronomical community. 



In addition to the possibility of detecting dark matter MACHOs, microlensing surveys 
also offer the interesting opportunity to study the populations of stellar lenses (which con- 
stitute at least a significant fraction of all events) and sources. Since most stars are located 
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in multiple systems, many microlensing events show deviations from the standard Paczyhski 
(point lens, point source) light curve. For example, some events show the effects of caustic 
crossings (Mao & Paczyhski 1991; Gould & Loeb 1992); this provides clear evidence of lens 
binarity and can be used to constrain the combination of the distances to the source and lens, 
the source-lens relative motion, and the physical extent of the source, as well as revealing 
the mass ratio of the lens system. A small additional fraction of microlensing events should 
show clear signatures of source binarity in multiple peaks and color shifts; this fraction was 
estimated to be 2-5% by Griest & Hu (1992). To date no convincing example of such an 
event has been reported in the literature; the best case is MACHO-96-BLG-4, but the light 
curve of this event can be equally well fit by a binary source or binary lens model, with the 
binary lens providing the most natural interpretation (Alcock et al. 2000). Several authors 
(Dominik 1998; Han & Jeong 1998; Di Stefano 2000) have proposed explanations for the 
lack of such clear binary source microlensing events. Essentially, a variety of different ef- 
fects including blending, unequal luminosity of binary components, and simple coincidence 
conspire to make most binary source microlensing light curves resemble single source, single 
lens events. A handful of other events, such as MACHO-96-LMC-2 Alcock et al. (2001), 
show signatures of binary orbital motion of the source; in this particular case, it is not clear 
whether both components of the binary contribute significantly to the total source flux. In 
any case the light curves of such events do not match the expectation of having multiple 
distinct, well-resolved peaks as described above. 

This work analyzes the light curve of OGLE-2003-BLG-095, which shows two distinct 
microlensing peaks. The light curve can be fit by a binary lens or a binary source model; 
however, the binary source model is statistically preferred. The light curve is presented in §2, 
§3 contains a description of the models and fitting procedure, and the results are analyzed 
in §4. 



2. OGLE data 

An /-band light curve of OGLE-2003-BLG-095 was obtained from the OGLE Inter- 
net archive 1 (cf. Udalski et al. 2003), and is shown in Figure 1. The source is located at 
(«J2000j ^J20oo) — (17:59:03.04, —27:11:19.8); a finding chart is also available from the OGLE 
web site. The data set consists of 109 observations from three observing seasons spanning a 
period of 2.2 yr, from August 2001 through 17 October 2003. The light curve from the 2001 
and 2002 observing seasons (seasons I and II hereafter) is relatively constant (to within ±3% 



1 http://sirius.ast rouw .edu.pl/~ogle/ 
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excluding the highest point and the lowest point), with a mean /-band magnitude of 15.58. 
However, due to the high stellar surface number density in the field, some blending seems 
likely. The light curve from the 2003 season (III) shows two very distinct smooth peaks, 
with maximum magnification factors of ~ 1.9 and m 1.7. For both peaks, the rise and fall 
in source intensity are well covered by observations. 

In order to simplify the fitting procedure, the light curve was converted to units of 
relative flux by adopting the mean intensity from seasons I and II as the base flux level. Thus 
the mean relative flux from seasons I and II is unity by construction. Even for relatively 
bright sources like the one in consideration, the photometric error estimates reported by the 
OGLE data reduction pipeline (cf. Udalski et al. 2002) are typically somewhat too small to 
account for the observed scatter in the light curves of sources that do not vary intrinsically 
(B. Paczyhski 2004, private communication). In the OGLE-2003-BLG-095 light curve from 
seasons I and II, the standard deviation s in flux is larger than the mean quoted error a by a 
factor s/a — 2.65, a ratio which nevertheless is significantly too large to be readily explained 
by underestimation of uncertainties. This is probably an indication that the source (or at 
least one component of it, if it is a blend) is variable at the level of at least a few percent. 
A careful inspection of the season II light curve indicates that variations are correlated from 
night to night. However, no periodic signal is present at a level of significance above 50% in 
a Lomb-Scargle normalized periodogram (cf. Press et al. 1992; Schwarzenberg-Czerny 1998) 
computed for seasons I and II. 

Such non-periodic variability substantially complicates the modeling process, causing 
greatest confusion for model effects that produce variations with similar amplitudes and 
time scales. See §4 for further discussion. To aid in interpreting fit results, the errors 
were scaled so that x 2 P er degree of freedom v (DOF), also known as the reduced x 2 an d 
denoted xli is unity for a constant fit to the light curve from seasons I and II. This required 
scaling all relative flux errors by a factor of 2.49. The scaling factor is similar to the ratio 
s/a quoted in the previous paragraph; the fact that these two quantities differ could either 
be attributed to a non-Gaussian error distribution (which seems likely given the apparent 
systematic variations) or to the limited number of data points. 

3. Modeling the light curve 

Since the intensity of the source during seasons I and II was taken to be constant, 
microlensing models were applied to the 50 data points from season III only. Two classes 
of models were considered: (1) models with a coplanar binary lens consisting of two point 
masses, in which the source was taken to be a small (but finite) disk of uniform surface 
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brightness ("binary lens" or "double lens" models); (2) models with a single point mass lens 
and two distinct point sources ("binary source" or "double source" models). The source 
was treated as finite in the former case due to the possibility of caustic crossings; although 
there is no direct evidence of such in the OGLE-2003-BLG-095 light curve, the possibility 
was considered in order not to restrict the exploration of parameter space (e.g., a case in 
which a brief caustic crossing event occurs in a gap between observations). In the latter case, 
such caution was not needed since the smooth, low-magnification light curve precludes any 
possibility of very small (close to zero) impact parameters. 

The light curve of OGLE-2003-BLG-095 suggests that the two binary components 
(whether of the lens or of the source) are separated by an angular distance of at least 
twice the angular Einstein radius 9e- For a typical value of 9e ~ 0.5 milli-arcseconds (mas) 
and a distance to the binary (lens or source) of at least a few kiloparsecs (kpc), which is very 
likely for an event in the direction of the Galactic bulge, the implied projected separation of 
binary components is several Astronomical Units (AU). Assuming the masses to be not too 
different from a Solar mass (M©), the minimum orbital period of the system is several years. 
Therefore, both classes of models were restricted to the static case; that is, any relative 
motion of binary components (lens or source) was neglected. 

3.1. Parallax effect 

For both classes of models, two parameters are needed to characterize the vector mi- 
crolens parallax 7r E (e.g., Gould 2000). One parameter can be chosen to be the parallax 
magnitude tte, whose inverse gives the projected size of the Einstein radius in the ob- 
server plane. A second parameter if) then gives the lens-source proper motion angle on the 
sky (measured counter-clockwise from North). Following the method recently proposed by 
Gould (2003), all parallax fitting is done in geocentric reference frames defined with respect 
to the position and velocity of Earth at times very close to one of the light curve peaks. 

3.2. Binary lens model 

Binary lens models are characterized by eight parameters, in addition to the two compo- 
nents of the microlens parallax. Four parameters describe the trajectory of the lens relative 
to the source: the time to of closest approach to the more massive binary component (corre- 
sponding roughly to one of the two event maxima), the Einstein time scale £e, the angle (f> 
between the lens line of centers and the lens trajectory with respect to the source position, 
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and the impact parameter uq to the more massive binary component. The remaining four 
parameters describe the lens and source properties: the binary mass ratio q and separation 
a, the fraction /y of the baseline intensity due to the source, and the ratio of the angular 
source radius to the angular Einstein radius, £* ; e = 0*/$e- All length scales (a, Uo) are in 
terms of the Einstein radius te of a lens with mass equal to the total mass of the binary. 

For events in which it has been possible to measure or estimate the ratio due to 
caustic crossings, typical values for microlensing toward the Galactic bulge are found to be 
on the order of 0.001-0.01 (Alcock et al. 2000; Jaroszyhski 2002). In the absence of caustic 
crossings and at low magnification, as in the case of OGLE-2003-BLG-095, binary lens model 
light curves are very insensitive to the precise value of £* 5 e provided it is smaller than a few 
tens of percent. Hence this parameter is poorly constrained by the available data, and its 
value is fixed at 0.01 in the subsequent analysis. 

For a binary lens the model flux can conveniently be found numerically for a given 
lens-source configuration by finding the (multiple) image locations for a set of points on 
the source boundary (e.g., Gould & Gaucherel 1997). This procedure gives 3-5 new sets of 
points defining the image boundaries, depending on whether the source is outside, crossing, 
or inside a caustic. Since the surface brightness of the source is taken to be constant, the 
magnification is given by the ratio of the sum of the (appropriately signed) areas enclosed 
by the image boundaries to the area of the source. 



The static binary source model is constructed from a linear combination of two point 
lens, point source models with Einstein time scales and vector parallaxes constrained to 
be the same for each source component. This model has seven parameters in addition to 
parallax: two for the fraction of the baseline intensity due to each source (fu,i, where i — 1,2), 
three for the trajectory (-u o ,£o?0 as m the previous case except that the reference point is 
defined by the closest approach of the lens to the source responsible for the first peak in the 
light curve), source separation a, and the Einstein time scale t&. 

The model flux m(t) at time t can be expressed as 



3.3. Binary source model 
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gives the total magnification for a point source, point lens system with separation u (in units 
of r E ). In this parametrization, the condition / bljl + f h \ 2 < 1 must clearly be met in order 
for the model to be physically meaningful. 

3.4. Fitting procedure 

The ^-minimization was performed using MINUIT, a function minimization package 
available as part of the CERN program library 2 . For each class of models, the fitting algo- 
rithm was initialized in several different realizations with different sets of initial parameter 
values. All such sets were chosen to produce an initial model light curve that was quali- 
tatively similar to the observed one (two well-separated smooth peaks with approximately 
the correct amplitude and width). For the binary lens model, four different initial source 
trajectories were chosen: two in which the source encountered the more massive lens first; 
two in which the source encountered the less massive lens first. For each of these pairs, one 
trajectory was chosen to pass in between the two lenses, and one was chosen not to pass in 
between the two lenses. The initial values of the other parameters were chosen appropriately 
to reproduce qualitatively the light curve shape for each of these cases. For the binary source 
model, an equivalent variety of initial parameter guesses were used, subject again to the re- 
quirement that the model light curve be qualitatively similar to the observations. After the 
best fit solutions were obtained for models in which the parallax was fixed at zero, full nine 
parameter fits (including parallax) were performed using the previous solutions as the initial 
guesses. 

3.5. Fitting results 

The best fit no-parallax model light curves and fit residuals are shown in Figure 2. 
Tables 1 and 2 contain the best fit results for the binary lens and binary source models, 
respectively. For the same number of DOF, the no-parallax double source model is preferred 
over the no-parallax double lens model by a highly significant margin of Ax 2 = 13.3, and 
furthermore this model satisfies the criterion that xl ~ 1 (with the caveat that the overall 
error scaling is rather uncertain). However, the lower panels of Figure 2 both show evidence 
of systematic residuals, indicating that unmodeled complexity is present in both fits. This is 
not at all surprising in light of the evidence for low-amplitude source variability mentioned 
in §2, and considering that parallax has been neglected. 



2 http:/ /cernlib. web. cern.ch/cernlib/index. html 
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Including the parallax effect adds two more free model parameters and results in Ax 2 = 
—25.7 and Ax 2 = —16.5 for the binary lens and binary source models, respectively. These 
large Ax 2 values formally indicate that the parallax effect is detected with a high level of 
statistical significance. The inclusion of parallax brings both models into the xl < 1 regime, 
for the adopted error scaling. It also reduces the difference in the statistical quality of the 
two model fits, though the double source model remains superior by A% 2 = 4.15 for the 
same number of DOF. 

It is worthwhile to note that although the overall error scaling is rather arbitrary as dis- 
cussed in §2 (and hence the precise value of x 2 is f ar from rigorous), Ax 2 is still a meaningful 
quantity. The low values of xl — 0.783 and xl — 0.682 respectively obtained for the best 
fit double lens and double source models indicate that the errors have probably been scaled 
by too large a factor, resulting in values of x 2 that are systematically low. In principle, the 
errors could be rescaled so that xl — 1 f° r the best fitting model, but this does not seem 
warranted, especially in view of the complications introduced by the likely source variability. 
Regardless, if the errors had been scaled by a smaller factor, A% 2 values would increase 
correspondingly. 

Additionally, there are a number of degeneracies in the problem that are important to 
consider. In the absence of parallax, that is, when the apparent lens-source relative trajectory 
is (modeled as being) strictly linear, there are two fully degenerate binary lens models and 
four fully degenerate binary source models. The degenerate models differ only by the sign of 
the impact parameter Uq and the orientation angle <fi. There are only two equivalent binary 
lens solutions because trajectories on which the two lens components pass on opposite sides 
of the source produce fundamentally different light curves (due to proximity to caustics) 
from trajectories on which the two lens components pass on the same side of the source. 
The same does not hold for the binary source models, and hence in this case there exist four 
degenerate trajectories. 

These degeneracies are somewhat alleviated when the parallax effect is included. For 
the binary lens model the degeneracy is broken; the two solutions separate by A% 2 = 2.79 
(only the preferred solution is given in Table 1). For the binary source model the degeneracy 
is largely broken as well; only the solution with the lowest x 2 is given in Table 2. The three 
solutions not included in Table 2 all fall within the range 1 < A% 2 < 2 of the preferred 
solution. Regardless, the differences among these near-degenerate solutions are largely im- 
material, since the best fit parameters (neglecting the sign of the impact parameter w an d 
the orientation angle 0) lie within or very near the 1-a error ranges quoted in Tables 1 and 2. 
For all the solutions, the parallax magnitudes 7Te are consistent within 1-a with the values 
from the best fit solutions. 
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There is a further possible degeneracy discovered recently by Gould (2003). This "jerk- 
parallax" degeneracy can produce multiple parallax solutions for a given sign of the impact 
parameter. In an effort to find any such solutions, a detailed search of the 7Te plane was 
performed for the double source model (which is more computationally tractable than the 
double lens model). The results are shown in Figure 3. No fully degenerate solution is 
present, although the contours are highly elongated in a direction perpendicular to the 
projected acceleration vector of Earth at the time of the first event maximum, when the 
geocentric reference frame is defined (as expected for this degeneracy). Numerical searches 
for other parameter combinations that produce similar trajectories were also performed (by 
solving eqs. 12-14 from Gould 2003) for both classes of models. No further degenerate 
solutions were discovered. 



4. Discussion 

An extremely relevant question is whether the parallax effect can be disentangled from 
the likely variability of the source (regardless of whether the variability is inherent or due 
to a hitherto unrecognized systematic effect). To address this issue, Figure 4 compares the 
deviations from constant flux in the season II light curve with the magnitude of the effect 
of parallax on the model light curves for season III. It is apparent that the modulations of 
the model light curves due to the parallax effect have similar amplitudes and variation time 
scales to the observed variability of the source, rendering the parallax measurements rather 
untrustworthy. 

Further support for the notion of "contaminated" parallax measurements comes from 
examining the implications of the best fit parallax-included models. There is a simple re- 
lationship between the projected Einstein radius te, the lens mass M, and the lens-source 
relative distance D re \ = (D~[ l — D^ 1 )^ 1 , where Di and D s are the distances from the observer 
to the lens and source; this relation is given by 

D = _£_1e_ (3) 

rel AGM 1 ; 

(Gould 2000). It is plausible to assume the source to be located in the Galactic bulge at a 
distance of approximately 8 kpc and to assume the lens to have a total mass of approximately 
one Solar mass. Substituting the measured values of 7Te from the best fit binary lens and 
binary source models, one obtains D m \ m 0.25 kpc and D rc \ 0.06 kpc, respectively. In both 
cases D~ c ] 3> D" 1 , indicating that D\ ~ A-ei- Even accounting for the substantial additional 
uncertainties introduced by assuming the lens mass and source distance, such a nearby lens 
seems unlikely. 
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Since the apparent variability of the source limits the degree to which the binary lens 
and binary source models can be statistically distinguished, a definitive interpretation of 
OGLE-2003-BLG-095 will probably require further observations. At the least, future OGLE 
observations of the source will allow better constraints on its variability. A simple (though not 
necessarily straightforward) test which might clearly confirm the binary source model could 
be performed through moderate resolution optical spectroscopy to see if the source exhibits 
double absorption lines or radial velocity variations. An alternative possibility would be to 
obtain deep, high resolution imaging of the field to see if the lens can be observed directly. 



The light curve of OGLE-2003-BLG-095 shows two well-separated smooth peaks that 
can be modeled as either microlensing by a binary lens or microlensing of a binary source. 
Both models are physically plausible and can possibly be tested by future observations of the 
source and/or lens. The binary lens model provides a more familiar explanation of the event 
since such phenomena have been observed many times in the past, while the binary source 
model is preferred on statistical grounds. Given the apparent low-amplitude variability of the 
source (or at least one component of it) however, statistical distinctions of the magnitude that 
separate the two models are somewhat suspect, especially considering the improbably large 
microlens parallaxes indicated by the best fit models. Despite these complications, OGLE- 
2003-BLG-095 remains a plausible candidate for microlensing of a well-resolved binary source 
- one of only two such observed to date. 

The author wishes to thank B. Paczyhski and T. Sumi for helpful comments and sug- 
gestions, M. Jaroszyhski for providing routines (originally developed by S. Mao) to calculate 
binary lens light curves, A. Gould for helpful discussion, and the OGLE collaboration for 
making microlensing data publicly available in real time. This research was supported in 
part by a NDSEG Fellowship award to the author. 
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Fig. 1. — /-band light curve of OGLE-2003-BLG-095 with nominal uncertainties from the 
OGLE data reduction pipeline. The observing seasons referred to in the text are labeled. 
The source is relatively constant during the first two seasons, with the deviations due to 
microlensing confined to the third season. 
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Fig. 2. — Panel (a): Best fit no-parallax binary lens (dotted curve) and binary source 
(dashed curve) models for the light curve of OGLE-2003-BLG-095. To the eye, there is little 
difference in the qualities of the two fits. Panel (b): Fit residuals in units of sigmas for 
the best fit binary lens model (not including parallax). Panel (c): Same as panel (b), but 
for the binary source model. The residuals in both of the lower panels display a somewhat 
systematic character. 
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Fig. 3. — x 2 contours in the 7r E plane for the double source model. 'X' marks the best fit 
solution; contours corresponding to A% 2 = 1, 4, 9, 16, 25, 36, 49, 64 are shown. The subscripts 
n, e refer to the North and East components of the parallax vector. The symbols || and _L 
refer to the projected direction of the Earth's acceleration vector at the time of the first event 
maximum (which is the reference point of the geocentric frame in which parallax fitting is 
done). Although the contours are elongated in the expected direction, no fully degenerate 
solution is found as in Gould (2003). 
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Fig. 4. — The curves do not represent fits to the data. Points with errorbars show devia- 
tions from constant relative flux from season II with nominal OGLE data reduction pipeline 
uncertainties. The dotted line is the binary lens parallax-included fit minus no-parallax fit 
from season III; the dashed line is the same, but for the binary source model. The qualitative 
similarities between the season II light curve and the season III "model-difference" curves are 
striking, especially considering that they refer to different observing seasons. The time scales 
for variation are similar, and notably, the amplitude of deviations from constant flux during 
season II is comparable to or greater than the differences between the parallax-included and 
no-parallax fits, even though no scaling has been applied to account for the fact that the 
source was brighter during season III (and thus relative flux deviations due to intrinsic source 
variability should be magnified as well). 
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Table 1. Parameters of best fit double lens model 

Quantity No-parallax Parallax included Units 

best fit value best fit value 

7T E 0.0 a 0.7 ±0.2 AIT 1 

i> 0.0 a —0.2+8:2 rad 

q 1.28 ±0.05 1.3 ±0.1 

a 2.90 ±0.07 3.1 ±0.1 r E 

0.048ig;88t 0.25l°;°? rad 

u -0.49 ±0.02 — 0.51±g;gi r E 

to 884.3 ±0.3 884.5 ±0.5 days b 

t E 39 ± 1 38 ± 2 days 

/.i o.92±g : gg o.95±8:8i 

£*,e 0.01 a 0.01 a 9 E 

X 2 57.76 32.10 

DOF 43 41 

X 2 U 1.34 0.783 



Note. - Parameter values for the no-parallax fit 
refer to a heliocentric frame; those for the parallax- 
included fit refer to a geocentric frame. Uncertainties 
are quoted at the A% 2 = 1 level. 

a fixed 



b HJD-2452000 
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Table 2. Parameters of best fit double source model 

Quantity No-parallax Parallax included Units 

best fit value best fit value 

7T E 0.0 a lAtoi AIT 1 

ip 0.0 a 0.07lg;^ rad 

/bi,i 0.28 ±0.03 0.34 ±0.04 

M 0.23 ±0.02 0.28 ±0.03 r E 

t 775.4 ±0.1 775.3 ±0.2 days b 

t E 44 ±2 41±| days 

/bi, 2 0.33 ±0.04 0.28±g;^ 

a 2.4 ±0.1 2.3+°;? r E 

-0.046 ±0.005 0.03lg;^ rad 

X 2 44.45 27.95 

DOF 43 41 

Xl 1.03 0.682 



Note. — Parameter values for the no-parallax fit refer 
to a heliocentric frame; those for the parallax-included 
fit refer to a geocentric frame defined by the position 
and velocity of Earth near the time to- Uncertainties are 
quoted at the A% 2 = 1 level. 

a fixed 



b HJD-2452000 



